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ABSTRACT 

We present measurements of the H i spin temperatmes (Tg) of the Cold Neutral 
Medium (CNM) towards radio sources that are closely aligned with stars for which 
published H2 ortho-para temperatures (Tqi) are available from UV observations. Our 
sample consists of 18 radio sources close to 16 nearby stars. The transverse separa- 
tion of the lines of sight of corresponding the UV and radio observations varies from 
0.1 to 12.0 pc at the distance of the star. The ultraviolet (UV) measurements do 
not have velocity information, so we use the velocities of low ionization species (e.g 
Na 1/ K 1/ C i) observed towards these same stars to make a plausible identification of 
the CNM corresponding to the H2 absorption. We then find that Toi and Ts match 
within observational uncertainties for lines-of-sight with H2 column density above 
10^^'® cm~^, but deviate from each other below this threshold. This is consistent with 
the expectation that in the CNM Tg tracks the kinetic temperature due to collisions 
and that Tqi is driven towards the kinetic temperature by proton exchange reactions. 

Key words: ISM: atoms - ISM: molecules - radio lines: ISM - ultraviolet: ISM. 



1 INTRODUCTION 

Physical conditions in the interstellar medium (ISM) have 
traditionally been studied using spectral lines from a variety 
of tracers including the 21 cm line and recombination lines of 
hydrogen in the radio regime, Lyman lines of H I, the Lyman 
and Werner bands of H2 and the atomic fine-structure lines 
such as C I in the UV as well as a host of rotational lines 
from molecules in the mm wavelength regime. For example 
the kinetic temperature of the gas can be determined using 
either the H I 21 cm line or the H2 UV lines, the pressure and 
cooling rates can be determined from the fine structure lines 
of C I and C 11* etc. Since many of these tracers co-exist in 
the diffuse ISM, multi-wavelength observations would allow 
one to cross check different observational techniques as well 
as to derive a more complete understanding of the physical 
state of the ISM. 

The gas temperature is a particularly well suited exam- 
ple of a parameter that can, in principle, be determined by 
observations at a variety of wavelengths. In the radio regime, 
the classical method consists of observing the H I 21-cm line 
in absorption towards a bright radio continuum source; this, 
in conjunction with observations of the emission spectrum 
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along a nearby line of sight allows one to measure the spin 
temperature (Tg) of the H I (see e.g. lKulkarni fc HeilesllOS^ . 
for details). While the H I spin temperature, strictly speak- 
ing, characterizes the population distribution between the 
two hyperfine levels of the hydrogen atom, it is often used 
as a proxy for the kinetic temperature of the gas. This is be- 
cause, in high density regions, Ts is expected to be tightly 
coupled to the kinetic temperature via collisions, while in 
low density regions, resonant scattering of Lyman-a pho- 
tons again may couple the spin temperature to the kinetic 
temperature (iField 1958). UV observations of the Lyman 
and Werner bands of H2 also allow one to determine the 
gas temperature. This is the so called "ortho-para" temper- 
ature (Toi). The ortho-para temperature characterizes the 
population distribution between the ortho and para forms 
of the H2 molecule, and like the spin temperature, it too 
is expected to be coupled to the kinetic temperature of the 
gas. This happens mainly via proton exchange collisions (e.g. 
LDakarno, Black & Woishoit 1973). In regions where H I and 
II2 co-exist, one might hence expect that the temperatures 
derived from radio (i.e. Ts) and UV observations (i.e. Tqi) 
should match. However, there have been very limited multi- 
wavelength studies of this sort. The main reason for this is 
probably that a given line of sight is rarely suitable for ob- 
servations at more than one wavelength. For example, UV 
observations are generally made toward bright nearby stars 
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which have no detectable radio emission. This rules out com- 
plimentary 21-cm studies along these lines of sight. Although 
a direct comparison of the two temperatures along a given 
line of sight is difficult it is still possible to compare the 
average properties of the two temperatures measured from 
different surveys. For exampl e, using a large sam ple of UV 
spectra towards nearby stars. ISavaee et al] (Il97'/i) found the 
mean value of Toi to be 77 ± 17 K, which is in good agree- 
ment with the mean value of the 21-cm Ts in CNM. This 
lends support to the argument that both these temperatures 
trace the kinetic temperature of the gas. 

The recent large scale high sensitivity and high angular 
resoluti on radio surveys lik e the NRAO VLA Sky Survey 
rNVSS. ICondon et alJll998h however do allow a way around 
the problem of finding lines of sight suitable for both UV and 
radio observations by identifying radio sources that happen 
to be close to the line of sight to UV bright stars. These 
radio sources are generally too faint for the classical sin- 
gle dish emission absorption studies, but are well suited for 
interferometers where the subtraction of the smooth back- 
ground emission is automatically achieved. In this work we 
present t he results of the Gi ant Metrewave Radio Telescope 
fGMRT. ISwarup et al.1ll991^ 21-cm H l observations toward 
18 radio sources. 



2 OBSERVATION AND DATA ANALYSIS 

Our sample (drawn from the NVSS catalog), consists of 18 
radio sources brighter than 100 mjy. For each source, we re- 
quire that there should be a sta r within 30' , a l ong th e line of 
sight to which Toi is available in lSavage et al.l (1197711 . Our 18 
radio sources correspond to 16 distinct stars, two stars hav- 
ing two corresponding radio sources. The details are summa- 
rized in Table 1. The columns in Table 1 are: (1) the name of 
the optical source and (2) the radio source, (3) the distances 
to the star, (4) the stellar co-ordinates, (5) the co-ordinates 
of the corre sponding radio sour ce, (6) the total H2 column 
density from lSavag e et aL* ('1977*), (7) the line of sight extinc- 
tion, E{B — V) from Savage et all il977ft. (8) the ortho-p ara 
temperature and its error bars ffrom ISavage et al]ll977l . as 
well as from our recalculations from their data, see below), 
(9) the H I 21-cm spin temperature (see sectionj^for details) 
and (10) the fractional difference between Toi and Ts. It is 
clear from Table 1 that lines of sight in our sample span 
wide range of E{B ~ V) and N(Il2). The transverse sepa- 
rations of the lines of sight of corresponding UV and radio 
observations at the distance of the stars ranges from 0.1 to 
12.0 pc. 

For the 6 sight lines with N(H2)^ 10^^ c m~^, no error 
bars f or the column densities were available in lSavage et alJ 
il977l) . For these cases, we obtained the column densities of 
H2 in J=0, and J=l rotational levels of ground vibrational 
level using Voigt profile fitted to the Copernicus archival 
data. Our e stimated column de nsities agree well with that 
reported bv lSavage et alJ lll977l) and we use our computed 
errors as indicative values while calculating the excitation 
temperatures for these systems. 

The GMRT radio observations were conducted between 
July 07-10 2002. For all sources the observing frequency was 
1419.4 MHz and the total bandwidth was 2 MHz with 128 
spectral channels (i.e. a velocity resolution of ~ 3.3 km s~^). 



The total on-source time was 3-4 hour on each source. Scans 
on standard calibrators were used for flux calibration, phase 
calibration and also to determine the bandpass shape. Data 
analysis was done using AIPS. After flagging out bad data, 
the flux density scale and instrumental phase were cali- 
brated. The continuum emission was then subtracted from 
the multi-channel visibility data set using the task UVSUB. 
Any residual continuum was then subtracted in the image 
plane by fltting a linear baseline to line-free regions using 
IMLIN. 

For the 21-cm emissio n spectra we took data fro m the 
Leiden-Dwingeloo survey jHartmann fc BurtonllT997l) . The 

angular resolution of this survey is ~ 36', i.e. larger than 
the separation between our target star and the correspond- 
ing radio source. The spectral resolution of the raw data is 
1.030 km s~^, however we Hanning smoothed these spectra 
to match the resolution of absorption spectra. 



3 CALCULATIONS 

3.1 Calculation of spin temperature: 

For a homogeneous cloud, the emission and absorption spec- 
tra uniquely yield the spin temperature, i.e. 



" 1.823 X 10i8/T(v)dv ^ ^ 

with N(H l) being determined from the (off source) emis- 
sion spectrum and t(v) being determined from the absorp- 
tion spectrum. In the real life situation where the gas is 
not homogeneous, but instead has density and temperature 
structure both along as well as transverse to the line of 
sight, the determination of "the" spin temperature from ra- 
dio observations is non trivial. As an example, application 
of Eq.Q to the observed spectrum produced by a set of op- 
tically thin multiple components along the line of sight will 
yield a column density weighted harmonic mean tempera- 
ture of the individual components. If the optical depths are 
large or there is structure both along and transverse to the 
line of sight, then in general there is no unique interpreta- 
tion of the data, althoug h, se veral app roaches to modeling 
have been attempted (e.g. jMebold 19721: iMebold et al.ll98a 
'Heiles & Troland 2003a b). We note that though there is no 
mathematically unique and physically robust procedure to 
interpret the spectra in this case, much of what we have 
learned about the neutral atomic medium has comes from 
the relatively simplistic assumptions underlying Eq. Q. In 
the analysis presented below, we continue in this tradition. 

There is in general absorption at several discrete ve- 
locity ranges arising from gas between us and the target 
star. As the absorption lines of H2 in J=0 and J=l lev- 
els are saturated, it is not possible to obtain the velocity 
information required to measure Toi for individual compo- 
nents; instead what is measured is the line of sight average 
value of Toi. It would hence be appropriate to also deal 
with the average Ts. Accordingly, we flrst independently 
decomposed the emission and absorption spectra into mul- 
tiple Gaussian components for the complete available ve- 
locity range using least-square minimization with minimum 
number of components. There is of course some measure 
of subjectivity in this decomposition. As a check on this. 
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Table 1. Details of our sample and results 



Star 


Radio Source 




Optical 




Radio 




log(NH,)° 


E(B-V)" 






Is 




A T/T 


HD 


NVSS 


(pc) 


(1, b) n 




(1, b) n 




_2 

cm 


mag 


K 






K 












22928 


J034227+474944 


82 


150.28, 


-05 


77 


150.20, 


-05 


78 


19.30 


0.01 


92 ± 


32 


98 


± 


13 


-0 


07 


± 


0.38 




J034629+474637 










150.76, 


-05 


41 










119 


± 


04 


-0 


25 


± 


0.35 


22951 


J034300+340634 


406 


158.92, 


-16 


70 


158.93, 


-16 


51 


20.46 


0.24 


63 ± 


14 


77 


± 


02 


-0 


20 


± 


0.22 


23408 


J034440+243622 


78 


166.17, 


-23 


51 


165.78, 


-23 


51 


19.75 


0.00 


89 ± 


40 


103 


± 


09 


-0 


14 


± 


0.45 


36486 


J053413-004408 


384 


203.86, 


-17 


74 


204.54, 


-17 


46 


14.68 


0.07 


1625+863'' 


219 


± 


41 


+1 


52 


± 


0.23 


36861 


J053450+100430 


532 


195.05, 


-12 


00 


194.89, 


-11 


98 


19.11 


0.12 


45 ± 


08 


92 


± 


03 


-0 


68 


± 


0.15 


37128 


J053550-012448 


409 


205.21, 


-17 


24 


205.36, 


-17 


42 


16.57 


0.08 


108 ± 


10'' 


108 


± 


05 


-0 


00 


± 


0.10 


38771 


J054409-091739 


520 


214.51, 


-18 


50 


213.75, 


-19 


14 


15.68 


0.07 


156 ± 


21'' 


72 


± 


01 


+0 


74 


± 


0.12 


40111 


J055703+261119 


1247 


183.97, 


+00 


84 


183.66, 


+00 


77 


19.74 


0.15 


117 ± 


52 


88 


± 


03 


+0 


28 


± 


0.44 


47839 


J064145+094704 


705 


202.94, 


+02 


20 


203.12, 


+02 


32 


15.55 


0.07 


1153+493'' 


54 


± 


02 


+1 


82 


± 


0.07 


57060 


J071741-241542 


1871 


237.82, 


-05 


37 


237.46, 


-05 


43 


15.78 


0.18 


82 ± 


06'' 


72 


± 


07 


+0 


12 


± 


0.13 


57061 


J071717-250453 


933 


238.18, 


-05 


54 


238.14, 


-05 


89 


15.47 


0.15 


513 ± 


46'' 


147 


± 


03 


+1 


11 


± 


0.06 


143275 


J160052-221214 


155 


350.10, 


+22 


49 


350.51, 


+22 


70 


19.41 


0.16 


56 ± 


12 


56 


± 


01 


+0 


00 


± 


0.22 




J160401-222341 










350.92, 


+22 


04 










48 


± 


00 


+0 


15 


± 


0.21 


147165 


J162400-261231 


142 


351.31, 


+17 


00 


351.28, 


+ 16 


11 


19.79 


0.38 


64 ± 


12 


53 


± 


03 


+0 


18 


± 


0.19 


149757 


J163631-105841 


138 


006.28, 


+23 


59 


005.82, 


+23 


47 


20.65 


0.32 


54 ± 


04 


80 


± 


07 


-0 


39 


± 


0.12 


209975 


J220320+624033 


1086 


104.87, 


+05 


39 


104.94, 


+05 


83 


20.08 


0.38 


77 ± 


21 


67 


± 


00 


+0 


14 


± 


0.26 


224572 


J000020+553908 


377 


115.55, 


-06 


36 


115.72, 


-06 


50 


20.23 


0.17 


82 ± 


23 


191 


± 


05 


-0 


80 


± 


0.24 



" ISavage et alHT977i) : ' Error in Tqi from our calculations 
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Figure 1. Profiles of the observed H I 21-cm emission 
and absorption towards 1=166°, b=— 23°.5 (corresponding to 
J034440+243622). The top panel is the emission spectrum taken 
from Leiden-Dwingeloo Survey (LDS). The bottom panel is the 
absorption spectrum taken using GMRT. Note that the broad 
component in emission is not detected in absorption. 

we note that although the decomposition was done inde- 
pendently for the emission and absorption spectra, for most 
cases we found components whose velocities match within 
the velocity resolution in the emission and absorption spec- 
tra. In some cases there is one or more wide components in 
the emission spectra which have no corresponding compo- 
nent in the absorption spectra (see Fig. 0. The traditional 
interpretation of this is that they arise from gas with too 
high a spin temperature to produce measurable absorption 
(Radhakrishnan ot al. 1972: .Heiles 2001: Hcilcs & Trolarfl 
l2003allbi) . i.e. the warm neutral medium (WNM), although 
of course it could also arise from a multiple weak narrow 
components or problems with the spectral baseline. Given 
the lack of velocity information in the H2 spectra, there is 



no unique way to match the H I absorption components 
found in this way with the H2 absorbing gas. We have in- 
stead used t he velocities of Na i, K i or C i absorption lines 
(taken from ("Jenkins, Jura & Loewenstein'lOSS'; 'Price et ^ 
2001; Wcltv, Hobbs & Kulkarni 1994; Wcltv & Hobbs 200^; 
Iwhite et al.ll200lf) 'l towards these same stars as a guide in 
this matching process. As the ionization potential of these 
species are less than that of H I and close to the energy re- 
quired to destroy H2, they are expected to coexist with H2 
in diffuse ISM. One should note however that (i) the com- 
bination of galactic rotation and velocity dispersion could 
cause distant H I to appear at the same velocity as nearby 
gas, (ii) there could be regions where these ions do not co- 
exist (i.e. the species are in a higher ionization state) but 
hydrogen is still in neutral state. Both of these mean that 
one cannot make a robust match between the H I and H2 ab- 
sorbing gas. We do find however that in most of the cases we 
can find matching components, and that these components 
by and large have |Vlsr| % 10 km s~^ (Given that the stars 
are relatively nearby |Vlsr| ^ 10 km s~^ is a plausible cutoff 
velocity for gas lying between us and the stars). We there- 
fore proceed by assuming that this matching is statistically 
correct. The average Ts was then calculated from the inte- 
grated emission and absorption spectra for all the matching 
CNM components. The calculation is done assuming that 
one is working in the limit of small optical depth, and also 
with no correction for absorption of the emission from one 
cloud by another cloud that happens to lie in front of it. This 
is for two reasons (i) as discussed above, there is no unique 
way to make this correction and (ii) as discussed in more 
detail below, we expect all such corrections to lie within the 
error bars of the Toi measurement. The errors in Ts were es- 
timated from the rms noise of the spectra and the FWHM of 
the Gaussian component for emission and absorption. This 
error should be taken only as an indicative value as it does 
not account for the error in fit and the error introduced 
through the assumption of small r. However, since the error 
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Table 2. Spatial variation of temperature: 
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Figure 2. AT/T as a function of different physical quantities. 
The horizontal lines in the figures mark AT/T = and 20% 
uncertainty. 



Source 


Position 


ds 


loeNl'Hl'l 


T 


HD 


1 (°) 


b (°) 


[pc] 


(cm~^) 


[K] 


HD22928 


150.28 


—05.77 






92 


radiol 


150.20 


-05.78 


0.12 


20.87 


98 


radio2 


150.76 


-05.41 


0.86 


21.14 


119 


HD23408 


166.17 


-23.51 


— 


— 


89 


radiol 


165.78 


—23.51 


0.48 


21.00 


103 


HD23480 


166.57 


-23.75 


0.60 




67 


HD36861 


195.05 


-12.00 


— 


20.78 


45 


radiol 


194.89 


-11.98 


1.45 


20.68 


92 


HD36822 


195.40 


-12.29 


4.16 


20.81 


63 


HD37128 


205.21 


-17.24 




20.45 


108 


radiol 


205.36 


-17.42 


1.67 


20.92 


108 


HD37742 


206.45 


-16.59 


9.65 


20.41 


101 


HD57060 


237.82 


-05.37 




20.70 


82 


radiol 


237.46 


-05.43 


12.00 


20.51 


72 


HD57061 


238.18 


-05.54 


12.95 


20.70 


513 


radio2 


238.14 


-05.89 


19.89 


20.24 


147 


HD143275 


350.10 


+22.49 




21.15 


56 


radiol 


350.51 


+22.70 


1.17 


21.09 


56 


radio2 


350.92 


+22.04 


2.38 


20.60 


48 



bars are in general small compared to the error bars in Toi , 
the indicative value should suffice. 



3.2 Calculation of ortho-para temperature: 

It is a standard procedure, in ISM studies, to use measured 
ortho-to-para ratio (OPR) of H2 to infer the kinetic temper- 
ature of the gas. If we assume most of the ortho and para H2 
reside in their respective ground levels then we can write, 

^ Nplij = 9 exp(-170.5/Toi) . (2) 

Toi derived from the above equation will either trace the ki- 
netic temperature of the gas (if OPR is controlled by proton 
or hydrogen exchange collisions) or formation temperature 
(if OPR is controlled by reactions in the grain surface) of 
H2 (f or details refer to ISternberg fc Neufeldlll99gl : IXakahasil 
120011) . We have used the observed column densities given in 
Table 1 and Eq.Q to estimate Toi. The associated error 
is computed using the standard error propagation. Since we 
can not resolve individual components, the derived Toi is 
the average temperature along the line of sight. 



and line of sight separation covered in our sample. How- 
ever, the points with AT/T > mainly come from the 
sightlines that have log N(H2)^ 15.8 or molecular fraction 
fH2(=2N(H2)/2N(H2)+N(H l)) ^ 10"''. This is expected as 
Solomon pumping dominates the excitation in this regime. 
The three points with AT/T < are from lines of sight that 
are optically thick in H2 with fna 5S 10~^ and may be region 
where H I and H2 are not co-existing. In panel (d) we plot A 
T/T as a function of maximum separation between optical 
and radio sightlines. We do not find any clear trend in this 
plot. We explore this issue further by collating in Table 2 the 
cases where we have more than 1 radio source or star within 
a separation of 1°. The columns in the table are (1) the 
name of the sources, (2) position, (3) separation of the lines 
of sight (ds) at the distance of the star, (4) N(H l) and (5) the 
temperature (Toi or Ts) measured along the line of sight. 
We find, among the cases where Toi — Ts, the match gener- 
ally improves with decreasing separation, though there are 
cases (e.g. HD36861) where the temperatures do not match 
even for a separation of 1.45 pc. One should note however 
that the separations are computed assuming the gas is at 
the distance of the star, while in reality the gas could lie 
anywhere between us and the star. 



4 RESULTS 



We define, the fractional difference between Ts and Toi by 
2(Toi - Ts) 



AT/T = 



(Toi + Ts 



(3) 



From Table 1 we can see that AT/T = within measure- 
ment uncertainties for nine out of 16 cases. For three sight- 
lines we find Ts > Toi and the other four sightlines have 
Ts < Toi . From the Fig. |5| it is also clear that the points 
with AT/T ~ spread over the whole range of E{B — V) 



5 DISCUSSION 

In this work we present the GMRT measurement of 21-cm 
spin temperature toward 18 radio sources that are close to 
16 bright starts for which UV observations are available. We 
find the Ts and Toi trace each other within the observational 
uncertainties when N(H2)^ 10^^ * cm~^ in 75 per cent of the 
cases. Toi is found to be higher than Ts for the sightlines 
with low H2 column density. 
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iHeiles fc Trolan have performed similar anal- 

ysis like us towards three directions. There is one star 
in our sample ( HP 22 951) that is in common with 
iHeiles fc Trolandl (Hoosi). For this sightline Toi = 63±14 
K and we find consistent Ts = 77±2 tow a,rd a radio source 
with a maximum separation of 1.35 pc. IHeiles fc Trolandl 
^2003a^ report Toi = 27±13 and 29±11 K toward NRAO 
140 and 3C93.1 that are separated by 14.6 and 9.5 pc re- 
spectively from the star. The spin temperature that they 
quote is that for a single CNM component along the line of 
sight. If we instead use the column density weighted average 
of all CNM com ponents in the veloc ity range probed by the 
C I absorption j Jenkins et alJlTgSSf) . the spin temperatures 
are 53 and 72 K, i.e. in agreement with Toi. 

A major concern is that we do not know if our tem- 
perature comparisons are for the same gas; we have only 
a plausible matching between the H I and H2 absorbing 
gas. As discussed above though, this matching is likely 
to be on the average correct. If this is so, then our find- 
ings mean that Ts and Toi do not vary much over dis- 
tances of a few parsecs. This is in apparent contrast to the 
well known findings that the H I 21-cm opacity varies on 
much smaller spatial scales (e.g . [Pigte r. Welch fc RomnevI 
Il97fil : ICrovisier. Dickey fc Kazes^ll^8R^ . The resolution mav 
be that the opacity fluctuations reflect not so much fine scale 
structure in the temperature as fine scale stru cture in the 
velocity and density fields (e.g. iBrogan et al.l .20fl5'l. It is 
also possible that our averaging over the several line of sight 
components decreases the effects of small scale variations. 

Another issue is related to the relatively large error bars 
in Tqi. Is it possible that all that our data is telling us is 
that the CNM has a characteristic temperature ~ 80 K, and 
that the general agreement between Toi and Ts merely re- 
fiects the fact that both the UV and radio observations are 
probing the CNM? Or does the temperature agreement ac- 
tually extend to individual lines of sight, as we have been 
asserting? A least-square linear fit using all 11 lines of sight 
for which Tqi and Ts agree (including both the radio sources 
close to HD22928 and HD143275) for Toi = rTs gives 
r= 0.981 ± 0.054. For all the lines of sight with H2 column 
density higher than our threshold, the least-square fit gives 
1= 0.778±0.080. As a quantitative check, we have calculated 
the Spearman correlation coefficient for increasingly larger 
subsamples. For a subsample with all 11 lines of sight for 
which Ts matches Toi, the correlation coefficient is max- 
imum (0.648). The corresponding significance is p < 0.05 
(from a two-tailed test). For a subsample with all 14 lines 
of sight for which II2 column density is more than 10^^'* 
cm~^, the correlation coefficient is 0.545. On the other hand 
when one includes lines of sight with N(H2) less than 10^® 
cm~^, the correlation coefficient goes down to 0.406 (for our 
complete sample). To summarize then, as per the Spearman 
rank coefficient test, at the better than 95% level, there is a 
one to one relation between Tqi and Ts for the 11 lines of 
sight with N(H2) higher than 10^^'* cm~^. Agreement be- 
tween Ts and Tqi for high column density gas would mean 
that the ortho-para equilibrium is mainly due to exchange 
collisions and that Toi does not reffect the formation tem- 
perature. The absence of relation between the two temper- 
atures when N(Il2) is optically thin is consistent with the 
slow rate of exchange collisions compared to the H2 destruc- 
tion rate and excess Solomon pumping from J=0 and J=l 



levels that make Toi deviate from the kinetic temperature 
and hence Ts. 
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